Codimension-two critical behavior in vacuum gravitational collapse 
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We consider the critical behavior at the threshold of black hole formation for the five dimensional 
vacuum Einstein equations satisfying the cohomogeneity-two triaxial Bianchi IX ansatz. Exploiting 
a discrete symmetry present in this model we predict the existence of a codimension-two attrac- 
tor. This prediction is confirmed numerically and the codimension-two attractor is identified as a 
discretely self-similar solution with two unstable modes. 



Introduction. Since the pioneering work of Choptuik 
on the collapse of self-gravitating scalar field [H, the na- 
ture of the boundary between dispersion and black hole 
formation in gravitational collapse has been a very ac- 
tive research area (see [2J for a review). One of the 
most intriguing aspects of these studies is the occur- 
rence of discretely self-similar critical solutions. Discrete 
self-similarity (DSS) means invariance under rescalings of 
space and time variables by a constant factor e A where 
A is a number, usually called the echoing period. Criti- 
cal solutions possessing this curious symmetry (with dif- 
ferent echoing periods) have been found for several self- 
gravitating matter models (massless scalar field [l| , Yang- 
Mills field [H, cr-model [1] and few others) and recently 
also in the vacuum gravitational collapse in higher di- 
mensions 

Our current understanding of DSS solutions is very 
limited in comparison with continuously self-similar 
(CSS) solutions. In the case of spherical symmetry the 
CSS ansatz leads to an ODE system which can be han- 
dled analytically and sometimes even rigorous proofs of 
existence are feasible. For example, the existence of 
a countable family of CSS solutions was proved in the 
Einstein-sigma model and the ground state of this 
family was identified as the critical solution (which was 
known previously from numerical studies of critical col- 
lapse). In contrast, the DSS ansatz leads to a 1 + 1 PDE 
eigenvalue problem which seems intractable analytically. 
Although this eigenvalue problem can be solved numer- 
ically.as was done by Gundlach for two models (scalar 
field [8| and Yang-Mills [9(), the numerical iterative pro- 
cedure requires a good initial seed in order to converge. 
Thus, Gundlach's method is efficient in validating and re- 
fining DSS solutions which are already known from direct 
numerical simulations but it is not useful in searching for 
new solutions. 

In this paper we consider the critical collapse for the 
five dimensional vacuum Einstein equations satisfying the 
cohomogeneity-two triaxial Bianchi IX ansatz and pro- 



vide heuristic arguments and numerical evidence for the 
the existence of a DSS solution with two unstable modes. 
This is a continuation of our studies in [f| where we have 
shown the existence of a critical DSS solution with one 
unstable mode and the associated type II critical behav- 
ior in this model. On the basis of our result it is tempting 
to conjecture that the critical DSS solution is a ground 
state of a countable family of DSS solutions with increas- 
ing number of unstable modes. 

Background. Our starting point is the cohomogeneity- 
two symmetry reduction of the Einstein equations in five 
dimensions which is based on the following ansatz intro- 
duced by us in || 



ds l = -Ae- 2S dt z +A- l dr l 
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where A, <5, B, and C are functions of time t and radius 
r. The angular part of (|TJ) is the 5J7(2)-invariant homo- 
geneous metric on the squashed 3-sphere with being 
standard left-invariant one-forms on SU (2) 
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-id9), a 3 = #-sin0 (ty. (2) 



where < 9 < tt, < (/) < 2tt, < ip < 4tt are the Euler 
angles. The squashing modes, B and C, play the role 
of dynamical degrees of freedom. This ansatz provides a 
simple 1 + 1 dimensional framework for investigating the 
dynamics of gravitational collapse in vacuum. In [5[ we 
made a simplifying assumption that B = C which means 
that the ansatz (JTJ) has an additional U(l) symmetry and 
only one dynamical degree of freedom (so called biaxial 
case). In this paper we drop this assumption and con- 
sider the full triaxial case with two dynamical degrees of 
freedom. 

Substituting the ansatz ([1]) into the vacuum Einstein 
equations in five dimensions and using the mass function 
m(t,r), defined by A = 1 — m/r 2 , we get the following 
system 
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where primes and overdots denote derivatives with respect to r and t, respectively. 



If B = C = 0, then the ansatz fl} reduces to the stan- 
dard spherically symmetric metric for which the Birkhoff 
theorem applies and the only solutions are Minkowski 
(5 = 0, m = 0) and Schwarzschild (<$ = 0, m = const > 
0). We showed in the biaxial case Q that these two 
static solutions play the role of attractors in the evolu- 
tion of generic regular initial data (small and large ones, 
respectively). We have verified that this property re- 
mains true in the triaxial case (see also [lQ(). Note that 
equations for small perturbations SB and SC around the 
Minkowski and Schwarzschild solutions decouple, hence 
linear stability results and the decay rates discussed in 
[B[ carry over without any changes. 

Heuristics. Below we focus on initial data on a bor- 
derline between dispersion and collapse. Our preliminary 
studies of the evolution of such data suggested that the 
relaxation of symmetry from biaxial to triaxial does not 
change the phenomenology of critical behavior observed 
in H , that is we have found our old biaxial DSS solution 
(hereafter referred to as the DSSi) as the critical solu- 
tion. In other words the biaxial symmetry seems to be 
recovered dynamically not only for generic initial data 
but also for the critical ones (see Fig. 1). The mecha- 
nism of this nonlinear synchronization phenomenon is an 
open problem which we hope to pursue elsewhere. Here 
we take it as a fact and use it as the starting point in 
further discussion. 

The second key element of our argument is the fact 
that the system (3) possesses a discrete Z 3 symmetry 
which corresponds to the freedom of permutations of co- 
efficients of one-forms 07. in the angular part of the metric 
((T|). These permutations are generated by the following 
transpositions 

T 12 : (B, C) - (C, B), T 23 : (B, C) - (B, -B-C), 
T 13 :(B,C)^(-B-C,C), (4) 

where the transposition swaps the coefficients of af 
and a 2 in JTJ). Biaxial configurations correspond to the 
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FIG. 1: Snapshots of the evolution of near-critical initial data 
of the form ([5]) for a — 1/2. In the last frame the squashing 
modes become synchronized and the solution coincides with 
the DS Si. 



fixed points of these transpositions: (B,B), {B, — B/2), 
and (B, -25). 

Thus, each biaxial solution, in particular the DSSi, ex_ 
ists in three different but geometrically equivalent copies. 
Let M. C rit denote the codimension-one critical surface in 
the phase space which separates dispersion from collapse 
and let Aii (i = 1,2,3) denote the basins of attraction 
of three symmetry-related copies of the DSSi solution. 
Since Mi lie in Merit and Merit is connected, there 
should exist codimension-two boundaries that separate 
different Mi from each other. It is natural to expect 
that these boundaries are given by the stable manifolds of 
three symmetry-related copies of a solution with two un- 
stable modes. In the next section we confirm this picture 
numerically and demonstrate that the codimension-two 
attractor is discretely self-similar. 
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Numerics. In order to find a codimension-two attrac- 
tor predicted above we consider the two-parameter family 
of time-symmetric initial data of the form 



B(0,r)=pf(r), C(0,r)=aB(0,r), 



(5) 



where f(r) is a smooth localized function satisfying reg- 
ularity conditions at r = 0. The results presented be- 
low were produced for the generalized gaussian /(r) = 
100r 2 exp(— 20(r — 0.1) 2 ). We pick some value of the 
parameter a and then, using bisection, we fine-tune the 
parameter p to the critical value p* (a) . We are interested 
in how the phenomenology of critical behavior depends 
on a. There are three distinguished values of a which cor- 
respond to three biaxial configurations: a = 1, a = —1/2, 
or a = —2. Since biaxiality is preserved during evolution, 
for these special initial data we obviously find our old 
DSSi as the critical solution (note that the correspond- 
ing critical values of p are related by the symmetry ([4]): 
p*(-l/2) = -2p*(l) andp*(-2) =p*(l) which provides 
a useful test of the accuracy of numerics). As we have 
already mentioned, numerical studies of critical collapse 
indicate that for other (generic) values of the parameter 
a the DSS\ solution, modulo the Z3 symmetry, also acts 
as an attractor. Let us denote by X^, X^', and X^ 
the three copies of the DSSi solution, corresponding to 
a = 1, a = — 1/2 and a — — 2, respectively. Below we 



focus our attention on X± and X^' ■ By continuity, the 
solution X± is the attractor for the values of a close to 
1 and the solution X% is the attractor for the values of 
a close to —1/2. Thus, in the interval —1/2 < a < 1 
there must exist at least one critical value a* such that 
the critical initial data (O corresponding to a* ± e (for 
a sufficiently small e) evolve to the different copies x[^ 
and X^~\ As usual, this critical value a* can be deter- 
mined by bisection We shall refer to initial data 
with parameters (a*,p*(a*)) as double critical. In Fig. 2 
we show the evolution of near double critical initial data. 
The key new feature which is apparent in Fig. 2 (in con- 
trast to Fig. 1) is the occurrence of another intermediate 
attractor around which the solution hangs for a while be- 
fore approaching the DSSi attractor. We find that the 
new attractor is discretely self-similar with the echoing 
period A 2 « 0.395 (see Figs. 3 and 4). We call this 
solution the DSS2 and denote by X^ 2 K 

The behavior seen if Fig. 2 has a natural explanation 
within the heuristic picture sketched above. According to 
this picture the DSS2 has two unstable modes, one tan- 
gential and one transversal to the critical surface. During 
the DSS2 intermediate regime the departure from X^ 
is well approximated by a linear combination of two un- 
stable modes 

6X^^ci(a,p)e^ >T fi(p 1 T)+c 2 (a,p)e^ >T f 2 (p,T), (6) 

where p = r / (T — t) and r = — ln(T— t) are the similarity 
variables (T is the blowup time of the DSS2), fi(p, t) are 
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FIG. 2: Snapshots of the evolution of near double criti- 
cal initial data (a = 0.1411..., p = 0.0952...). The new 
(codimension-two) intermediate attractor (DSS2) is seen in 
the second frame. The departure from this attractor (proceed- 
ing from the origin) and the approach to the old (codimension- 
one) intermediate attractor (X± copy of the DSSi) in seen 
in the third frame. In the last frame the DSSi regime is well 
developed near the origin while the DSS2 regime is still vis- 
ible for large r. To produce this figure we had to fine-tune p 
to the critical value with a relative accuracy of 10" 26 (using 
quadruple precision arithmetics). 
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FIG. 3: Evidence for discrete self-similarity of the 
codimension-two attractor. We plot the mode B at some 
central proper time To during the intermediate regime shown 
in the second frame in Fig. 2 and superimpose the next four 
echoes which subsequently develop. 

the eigenfunctions (periodic in r with the period A2) and 

(2) 

> are the respective eigenvalues. The coefficients 
Cj are the only vestige of initial data. Given a close to 
a*, for the perfect fine-tuning of p to its critical value, 
the transversal unstable mode (let it be /1) is completely 
suppressed, that is ci(<z,p*(a)) = 0. In practice, we can 
fine-tune p with very high precision so that ci(a,p(a)) is 
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FIG. 4: The close-up of the DSS2 solution. It is evident that 
this solution is not biaxial. Note that the mode C oscillates 
with twice the frequency of the mode B. 

almost zero. For such data C\ <C C2, hence after several 
cycles the solution leaves the DSS2 along the tangen- 
tial unstable direction and then, after a short transient 
period, approaches one of the copies X?> or (de- 
pending on the sign of C2) of the DSSi. Afterwards the 
solution stays near DSSi for some time and eventually 
moves away from the critical surface M. C rit towards col- 
lapse or dispersion. 

The existence of the DSS2 solution has an interest- 
ing consequence for the black-hole mass scaling law, 
AIbh ~ (p — P*) 7 , associated with the type II critical 
collapse. Consider a near critical value a and monitor 
the black hole mass over a wide range of supercritical 
values p > p*{a). As we have discussed above, for p very 
close to the threshold the solution moves away from the 
critical surface along the unstable mode of the DSSi. 
However, as p — p* increases, we observe a competition 
between the transversal and the tangential modes of the 
DSS2- Far from the threshold the transversal mode be- 
comes dominant (c\ 3> c%) which means that the solution 
moves away from the critical surface along the transver- 
sal unstable direction of the DSS2 and never approaches 
the DSSi. Accordingto the well-known argument from 
dimensional analysis [2J the scaling exponent 7 is related 
to the eigenvalue A of the unstable mode along which 
the solution is ejected towards collapse: 7 = 2/ A (the 
factor of 2 comes from the fact that in five dimensions 
mass has dimension length 2 ). Thus, for near double crit- 
ical initial data the mass scaling law is expected to ex- 
hibit the crossover from the exponent 71 = 2/A^ 1 ^ (where 
A^ 1 ) is the eigenvalue of the single unstable mode of the 

(2) 

DSSi) for small p — p* to 72 = 2/ A]/ for relatively large 
p — p* (see Fig. 5 for the numerical confirmation of this 

prediction). Using this relationship and the numerically 

(2) 

determined 72 we estimate that X\ ' m 6.67. The second 

(2) 

eigenvalue A 2 associated with the tangential unstable 



mode of the DSS2 can be obtained from the relationship 

(2) 

A 2 A2<5n = — Stt, where Sn and Sir are the changes in the 

echo number n and n = In \p — p* \ . Fitting this formula 

(2) 

to numerical data we obtain A 2 ~ 6.03. 



-4 



-8 



- 


1 1 1 


1 1 1 














7! ~ 0.328 






y 2 ~ 0.300 


y 

-.? 










.■■jf 


data 






linear fit 1 






linear fit 2 















-70 -60 -50 -40 -30 -20 -10 

ln(p-p* ) 



FIG. 5: The black-hole mass scaling law for near double crit- 
ical initial data. The linear fit yields different slopes for small 
and large distances from the threshold. 

We remark that the heuristic argument for the exis- 
tence of a codimension-two attractor presented in this 
paper can be applied repeatedly to argue for the exis- 
tence of higher codimension attractors [12| , however the 
numerical search for these solutions via multi-parameter 
fine-tuning would be extremely difficult. 
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